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DEMONSTRATION OF THE DIFFERENTIAL EQUATIONS 
EMPLOYED BY DELAUNAY IN THE LUNAR THEORY. 



BY G. W. HILL. 

The method of treating the lunar theory adopted by Delaunay is so ele- 
gant that it cannot fail to become in the future the classic method of treating 
all the problems of celestial mechanics. The canonical system of equations 
employed by Delaunay is not demonstrated by him in his work, but he 
refers to a memoir of Binet inserted in the Journal de V Ecole Polyteohnique, 
Cahier XXVIII. Among the innumerable sets of canonical elements it 
does not appear that a better can be selected. These equations can be estab- 
lished in a very elegant manner by using the properties of Lagrange's and 
Poisson's quantities {a, b) and [a, 6]. But a demonstration, founded on more 
direct and elementary considerations, is, on some accounts, to be preferred. 

Let a denote the mean distance, e the eccentricity, i the inclination of the 
orbit to a fixed plane, I the mean anomaly, g the angular distance of the low- 
er apsis from the ascending node, h the longitude of the ascending node 
measured from a fixed line in the fixed plane, // the sum of the masses of 
the bodies whose relative motion is considered, and R the ordinary pertur- 

bative function augmented by the term jrC_. Then if we put L = '^(jia), 

G- = Vlj"ffl(l — ^)\ H = Vljua(l — e*)]cos^, Delaunay 's equations are 
dL^dR dG _ dR 

dt dl' dt dg' 

dt~ dL' dt dG' 

In terms of rectangular coordinates 
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In this expression for «, y, z ought to be substituted their values deduced 
from the formulas of elliptic motion, and expressed in terms of i, G, H, 

I, g, h. It should be noted that the term -^^ == jp, of the zero order with 

respect to the disturbing forces, has been added to R only to preserve in the 
equations the canonical form : it is only by amplifying the signification of 
the word that I can be called an element, as it is not constant in elliptic mo- 
tion but augments proportionally to the time and — = n = -t-. It is cho- 
sen as a variable in preference to the element attached to it by addition 
simply to prevent t from appearing in the derivatives of R outside of the 
functional signs sine and cosine. 
The equations 

d? "^ r3~ dx' W^ r^ dy' de "^ r^ dz' 
are well known; here, however, R does not contain the term-c^. By mul- 
tiplying them severally by dx, dy, dz adding and integrating is obtained 
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Dividing both members of this equation by — n = — ^{fta~^), the left 
member is seen to be the differential of \/{fia) = L. Consequently 

dL _ dR 
dt dl ' 

Denoting the true anomaly by v, the orthogonal projection of the radius 
vector on the line of nodes is r cos (v+g), on a line perpendicular to it and 
in the plane of the orbit r sin (v+g). And the latter, projected on the plane 
of reference, is r8in(v-\-g)cosi, and on a line perpendicular to this plane 
r sin {v+g) sin i. If the two projections lying in the plane of reference are 
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again each projected on the axis of x, their sum will be the value of the co- 
ordinate X, and the sum of their projections on the axis of y, the value of 
the coordinate y. Hence 

X = rcos{v -\- g) cosh — r sin {v + 9) cos i sin h, 
y = rcos{v + g)sia.h + >■ sin {v -\- g) cosicos j^, 
z = rsin{v + g) sin i, 
or substituting for i its value in terms of G and H, 

X = rcos (v + S')cos h — -=f-rsin(v -f- g)smh, 

TJ" 

2^ = 7- cos (v + 5^) sin A + y^r sin (v + g) cos A, 

Or 

z = -^-i — — i r sm [v + g). 

As r and v are functions of L, G and Z only, the preceding equations show 
the manner in which H, g and h are involved in R. 

ff denotes double the areal velocity projected on the plane xy, or 

xdy — ydx jr 

dt 

Co„.e,uenay f=»f-!'§- 

But the foregoing values of x, y, z show that we have 

dx dy dz ^ 

dh~ y' dh~''' dh~ ' 

and thus ^ == ^^ + ^^ + i^^ = ^. 
dt dx dh dy dh dz dh dh' 

G denotes double the areal velocity, and evidently, if for the moment we 
suppose X and y to be drawn in the plane of the orbit, the axis of x towards 
the lower apsis, 

dt dy dx dv ' 

where, in the last R, for x, y and z must be substituted their values given 
above in terms of r, v, G, H, g, h. Now as the only way, in which g is 
involved in these values, is by addition to v, it follows that 

dR ^ dR , 
dv dg ' 

and this equation is not affected when, for r and v in R, are substituted their 
values in terms of L, G and I. Consequently 

dG _ dR 
dt dg ' 
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In the elliptic theory 

at 

Whence we deduce 

d[i/(ffl — g^)cos;t] ^ a; — — 2— 

dt dz dx' 

d\y{G' — H^)Bmh'\ ^ ydB_^dB 

dt ^ dz dy' 

Eliminating ' ^ . froni th^^se equations, we obtain 

dh _ 2 sin A dB z cos A dB a; sin A — y cos h dB 

'dt ~ VW^^^SPj ~d^ ^{G'—H') dy V ( Gr^—W) dz ' 
Comparing the coefficients of the three derivatives of B in the right 
member of this equation with the values of a;, y and z in terms of r, v, G, 
JET, £r, A, we recognize that they are severally equivalent to the negative of 
the partial derivatives of these quantities with respect to fl". So that 

dh i dB dx . dB dy . dB dz \ dB 

^~ [dx dH^ dy dH~^ dz dH) dH' 

It is a well known principle in the theory of varying elements, that if 
we differentiate any function, which is a function of the coordinates and t 
only, but expressed in terms of t and the elements, with respect to t only in- 
as much as it is explicitly involved, we obtain the correct value. Hence if 
the differentiation is performed on the supposition that the elements arc 
alone variable, the result should be zero. Applying this to the function r, 

we get 

dr^dL drdG,dr(dI__\ ^q 
dL dt ^ dG dt ^ dl \di I ' 



or 

or again 



dr dB ^ dB , dr^(di_ „\ __ a 

dLdf^dG da^ dl\dt ) ' 



dr(dBd>i,dBdv\ , dr d^ di^/dl_ _ ^\ _ q 
dL\dr dl'^ dvdll "^ dG dv "^ dl \dt I 



dL\dr dl 

Whence we derive 

dl _ drdB_ fdry^ fdr dv_ , dr'ldB 

-di~^~dLd? \dl) l.dLdl'^ dGJdv' 

From the expression for r we can eliminate I and introduce v in its place by 

means of the expression for v in terms of L, G and I; the result is the 

well-known equation 
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r- <-^-^) == ^! 

1+ecosv r-i^ViL^—O') 1* 

And we have 

dr _ dr^ dr _ / dr\ ,dr do 

dt^d^df' dL~\dL)'^dkdV 

the parentheses denoting the derivative with respect to L only insomuch as 

it enters the preceding equation for r. By making these substitutions, the 

coefficient of -j- in the expression for -j- becomes 

dv fdr\~^T~(d/r\dv idr'~\ 

dZ WJ \_\dL}W^dGfS 

From the preceding equation for r, we derive 

(dr\ fj.r' COS V 

\dL) ~ ''IN' ' 
also the following is a well known equation in the elliptic theory 

d^^Q_ 
dl rw-2' 

For obtaining the value of -r^, u being the eccentric anomaly, we have 

the equations r = a{l — ecos w), 1 = u — esinw. 
Their differentials give 

dr , . du 

-i- = — a cos M + ae sin M -y-, 
de de 

= (1 — e cos u)du — sin u.de. 

TTTi dr cosu — e 

Whence j- = — « :; = — a cos v. 

de 1 — ecosM 

And ,- ViL'-G') ^-_G 

dr dr de _ O cos v 

dG- ~ de dG ~ [xe 
By substituting the values, it is found that 

/dr\dv^ , dr r. 

[dljW^ dG~ 

T dl dr dR do dB 

In consequence -^=zn == -, y= t- 

dt dL dr dL dv 

_^_dR 
~ dL' 
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As i2 is a function of the coordinates and the time only, we can treat it 
as we have done r. Then 

dB dL ,dB(dl_ \ _._ dR dG .dRdg dB dH , dBdh_^ 
dL dt '^ dl \dt ^} '^ do dt '^ dg dt ^ dH dt '^ dh dt 
On substituting in this the values of the differentials of the elements which 
have already been determined, it is seen that all the terms but two mutual- 

dB 

ly cancel each other. And on dividing the result by -5—, we get 

dg dB 

dt~ dG' 

By adding to B the term ^^ = ^, its partial derivative with respect to 

„2 

L is augmented by the term — 73= — ^t but all the other derivatives are 

unchanged. In consequence of this addition the value of the differential of 

, , dl dB 

I becomes -r? = rr-- 

dt dL 

An objection may be made against the preceding method of obtaining the 
differentials of I and g, in that the quantities -^j and -=— , which both peri- 
odically vanish, have been employed as divisors. But this objection has 
force only when it is admitted that the differentials of I and g or the corres- 
ponding derivatives of B may be discontinuous. For having proved the 
truth of the equations for all times, except when the divisors, just mentioned, 
vanish, it follows, that if both members are continuous, the equations must 

dr dB 

still hold even for the moments of time when -7^ = or -,— = 0. 

dl dg 



BECENT BE8ULTS IN THE STUD Y OF LINKAGES. 



BY PEOF. W. W. JOHNSON, ANNAPOLIS, MD. 

{Contimied from page 46.) 

The so called kite-shaped quadrilateral, that which has two pairs of equal 
adjacent sides, may be employsd as a cell in which the rays are always at 
right angles. Let ABCD be such a quadrilateral the length of its unequal 
sides being a and b. Take any point on AB for the fulcrum and draw 
OQ and OP parallel to the diagonals. It is evident that however the 



